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Abstract. In this work, we investigate in detail the capabilities of present (H.E.S.S., 
MAGIC, VERITAS) and planned (CTA) ground-based Cherenkov telescope systems 
to detect angular anisotropies in the diffuse gamma-ray background. We first study 
the impact of instrumental characteristics (effective area, field of view, angular reso- 
lution, and background rejection efficiency) to the ability to detect anisotropies. In 
addition, we compare different observation strategies, i.e., whether a single deep obser- 
vation or a splitting over multiple shallow fields is preferred. Secondly, the sensitivity 
to anisotropies generated by self-annihilating dark matter is investigated for different 
common dark matter models. With planned configurations of CTA, we find that a rel- 
ative contribution of ~ 10% from dark matter annihilation to the diffuse gamma-ray 
background can be detected, together with the sensitivity to the self-annihilation cross 
section {av) = 3 x 10~^^ cm^s~^ expected from thermal dark matter freeze-out. We 
also stress the importance of constraining anisotropies from unresolved astrophysical 
sources already with the current generation of instruments, as a novel and complemen- 
tary method to constrain the properties of TeV sources. 
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1 Introduction 

The study of gamma-ray anisotropies [1, 2] is recently providing new and complementary 
insights into the nature of gamma-ray sources and the diffuse gamma-ray background 
(DGB) [3, 4]. The pattern of anisotropies has been studied with different techniques, 
mainly through its Angular Power Spectrum (APS), as in [1, 2]. Likewise, the study of 
the 1-point Probability Distribution Function (PDF) [5] and the cross-correlation with 
galaxy catalogues [6] provides further complementary information. 

It has been argued that Dark Matter (DM) self-annihilation or decay can leave a 
peculiar imprint into the anisotropy pattern and spectrum of the DGB [7-18]. In fact, 
while the emissivity of ordinary astrophysical sources scales with their inner densities 
~ Q (modulo a source-class dependent bias factor), the emissivity of self-annihilating 
DM scales with its density squared Owing to this difference, self-annihilating DM 

could leave a characteristic imprint in the anisotropy power spectrum of the DGB. This 
simple picture, however, can be spoilt by the presence of unresolved point sources which 
produce a Poisson-like APS, more closely resembling the DM one. Nonetheless, even if 
the astrophysical Poissonian term dominates the intrinsic clustering APS and, thus, the 
astrophysical and dark matter APS look indeed similar, it is still possible to separate the 
two contributions by looking at the energy dependence of the APS, i.e. looking at the 
anisotropy energy spectrum [13]. The different hypotheses about the origin of the DGB 
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can thus be tested by measuring the angular power spectrum together with its energy 
dependence. 

At the same time, theoretical predictions for the anisotropics generated by different 
gamma-ray source populations (e.g., blazars and galaxy clusters [19], millisecond pul- 
sars [20], star-forming galaxies [21]) are becoming available, contributing to clarify and 
extend our knowledge of this approach and its potential. 

At present, gamma-ray radiation is detected mainly through two techniques. Ob- 
servatories in space, such as the Fermi Large Area Telescope (Fermi-LAT) [22], detect 
gamma-rays through pair conversion in the detector itself. Fermi-LAT observes gamma- 
ray photons in the energy range from ~ 20MeV up to a few hundred GeV, with an 
effective area close to 1 m^ and a field of view (fov) of ~ 1 sr. Fermi-LAT routinely op- 
erates in sky surveying mode, continuously mapping the entire sky. On the other hand, 
ground-based Cherenkov telescopes like the current experiments H.E.S.S. [23], MAGIC 
[24], and VERITAS [25], or the planned Cherenkov telescope array (CTA) [26, 27] ob- 
serve the Cherenkov light produced by air showers initiated by gamma-rays in the upper 
atmosphere. The effective collection area of ground-based instruments is typically of the 
order of 10^-10^ m^. The energy range of the current experiments covers about 60 GeV 
to 100 TeV, but future realizations of this concept will lower the threshold to 10 GeV 
or even 5 GeV [28]. Contrary to Fermi-LAT, ground-based instruments cover only a 
relatively small field of view of typically a few msr, so that all-sky scans are not feasible. 
Observations are possible only during the night and in conditions of clear sky reducing 
the duty-cycle to about 10%. Furthermore, the trigger rate is dominated by a large back- 
ground of hadronic showers. Specific techniques are employed to reduce this background 
substantially. However, even after sufficient gamma-hadron separation, air showers in- 
duced by cosmic-ray electrons can still substantially contribute to the background at a 
few hundred GeV [29-32], hardly separable from photon-induced air showers. Despite 
the above mentioned difficulties, interesting scales for the investigation of anisotropics 
are typically very small (less than 1°), so that the small field of view is not a serious 
limit for their study. The angular resolution of Cherenkov telescopes is typically better 
than 0.1°, so that multipoles in the range 100-1000 can be easily resolved. At the same 
time, the background is expected to be isotropic at small scales and thus also not a 
fundamental obstacle. 

In this paper, we investigate in detail the capabilities of ground-based Cherenkov 
telescopes to measure gamma-ray angular anisotropics. In Section 2, we introduce a 
simplified Monte-Carlo approach, that we use to study the impact of instrumental char- 
acteristics (fov, angular resolution, and background rejection power) to the detection 
sensitivity for anisotropics. Since it would be very effective to analyze all the data col- 
lected from several observations taken for other scientific goals, we then investigate the 
influence of the observation strategy - whether it would be advantageous or disadvan- 
tageous to split the data into several fov. In Section 3 we list the instrumental setups 
that we use as benchmarks and we provide estimates of the CR background expected for 
these configurations. In Section 4, we use a more realistic simulation setup and study 
the sensitivity to anisotropics from DM self-annihilation using different DM models. We 
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provide final comments and discussion in Section 5. 



2 Optimizing the observation strategy 
2.1 Simplified setup 

As a measure of the anisotropy we will use the angular power spectrum (APS) of 
fluctuations throughout this paper. Given a map !{'&, ip) on the sphere, the fluctu- 
ation map is defined as 5{'d,{p) = I{-d,ip)/I — 1, where I is the mean value of /. 
The fluctuation map 6{'&,ip) is then decomposed into spherical harmonics Y^{'d,ip) as 
6{i!},ip) = Ylem'^irnY^^i^, where are the coefficients of the spherical harmonic 
decomposition and m = —I, From the aim, the APS is defined as 

Ce = {\aem\^), (2.1) 

where (...) indicates the statistical ensemble average. Then, the quantity 

C y h!Z^ (2.2) 

m 

is an unbiased estimator of the true power spectrum C^, i.e. (C^) = C^. 

Note that instead of the fluctuation dimensionless APS, the dimension-full APS of 
anisotropics (the APS of the map /(i?, ip) itself) can be used. This is particularly useful 
when analyzing real data (see [1]), but it will not be considered in the following. 

In order to simulate a list of events containing anisotropics, we generate sky maps 
with a given spectrum. Ci can be interpreted as the width of the aim distribution 
over m for flxed i. Assuming Gaussian fluctuations, the aim for a fixed value of I can be 
randomly chosen from a Gaussian distribution with mean and width y/Ci. The phase 
is chosen equally distributed between an 27r with the condition aim = (^*-im^ to ensure 
fluctuation maps with real values. For a given spectrum, we generate twelve realizations 
of the aim and thus twelve independent fluctuation maps of 5{'&,Lp). Moreover, we 
simulate maps with five different APS following a power law with slopes s = 0.5, 1.0, 
1.5, 2.0, 2.5, i.e. i{£ + l)Ci ~ The angular resolution of the simulated sky maps 
is chosen to be 0.002°, corresponding to a maximum resolvable multipole £ = 9 x 10^. 
The maps are normalized as (5'('t?, ^p) = (5(i?, p>) — (5mm) / ((^max — Smm), in order to obtain 
a distribution between and 1 to simulate events. Therefore, the map 6' acts as an 
effective intensity map /'(i?,(/3) = 5' (??,(/?). Note that the anisotropics of this effective 
intensity map are independent of the original 5{'d,ip) normalization and always in the 
range — 1 < 61' /I' < 1, thus giving large fiuctuations of the order of 100%. The setup has 
been optimized for the purpose of investigating the effects of instrumental characteristics 
on the APS. We will, however, extend the setup in Section 4, where we use a more realistic 
simulation of anisotropics, allowing to choose the desired anisotropy level. 

Based upon the template intensity maps !'{'&, (p), the simulation of events requires 
a choice of three parameters: (Tfov, the half-width of the camera acceptance, fipsf, the 
width of the point-spread function (PSF), and the signal fraction /sig, the fraction of 
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signal events with respect to the sum of signal and background events. Event lists 
will contain a signal component (anisotropic) and a background component (isotropic 
by definition). Both, the camera acceptance as well as the PSF are assumed to follow 
Gaussian distributions. 

Each event of the list is handled in three subsequent steps: The celestial position 
is chosen randomly according to the fov. Comparing a uniform deviate with /sig, it is 
then decided whether the event is handled as a signal or a background event. If the 
event belongs to the background, the event is just retained. In the signal case, instead, 
a variate z is generated following a normalized uniform distribution. If z is smaller than 
/'(t?, f) at the event's position, the event is kept, and otherwise is rejected. Finally, each 
event is randomly displaced from its original direction according to the PSF, to realize 
the convolution of the map with the PSF. If not particularly specified, an event list 
will contain A'event = 10^ entries. We remark that this number is unrealistically high, 
even for CTA, and is just used to isolate and emphasize instrumental effects. Realistic 
numbers of signal and background events are provided in the next section. 

We use the HEALPix software and pixelization scheme [34] to create and analyze 
count maps (with A'pix pixels), and to extract the APS. Before the analysis, the count 
map is cast into a fluctuation map 



where Xi denotes the number of events in pixel i, and bi{'d, if) equals 1 inside pixel i and 
outside. In this way, 5"{'d, ip) is normalized such that J dO, (1+6" {"d, ip)) = 4tt. Strictly 
speaking, 6"{'&,ip) does not represent a fluctuation map, a factor of fsky = ^/oii/47r 
being required for the correct normalization. However, the above normalization ensures 
that the derived APS is already corrected for the incomplete sky coverage effect, at least 
at high multipoles [1]. Also note that we are directly analyzing the fluctuation map 
build out of the raw count map, while, in principle, the fluctuation map build from the 
reconstructed flux map (the count map divided by the exposure of the experiment or 
the fov in this case) can be analyzed. While the latter would be preferable for a real 
data analysis, the former is sufficiently accurate for a sensitivity study. It implies the 
extracted APS to be a convolution of the experiment's windowing function (the APS of 
the camera acceptance) with the true signal. Moreover, for the large I considered below 
the windowing effect is marginal. 

The statistical uncertainty of the APS is derived from the twelve simulations per- 
formed. In the following plots, the mean and the RMS as estimate of the standard 
deviation are shown as uncertainty bands. 

2.2 Influence of the detector configuration 

In Fig. 1, examples of recovered APS derived from the maps simulated with a given input 
APS slope are shown, for a varying number of total events. No background is assumed 
and all events are of signal type. For small i, the windowing function dominates and the 
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Figure 1: Angular power spectra recovered from input event maps simulated -with 
different APS slopes s from 0.5 to 2.5 in steps of As = 0.5 (from red to blue). In this 
figure, PSF effects are neglected and crfov = 2.5°. No background is included in the 
simulations. Subfigures represent the cases of (a) 10'^ events, (b) 10^ events, and (c) 10^ 
events. The effect of an increasing event number on the Poissonian noise is depicted in 
(d), -which shows the APS of background-only simulations. 



spectrum is distorted. At large i and for a sufficiently large number of events, the signal 
starts to dominate, and different signal slopes can be discriminated. If the number of 
events is comparably small, however, the random Poissonian noise from the finite number 
of events will dominate the signal. The Poissonian noise is given by Cn = 4:Tt/N for full 
sky coverage (see appendix B in [9] and section IV. A in [7]), where N = Ngig + N^^g is 
the total number of events. This is in good agreement with the noise estimates from 
the simulations of background only events (and partial sky-coverage corrected APS, see 
previous section), as shown in the lower right panel of Fig. 1. 

The effect of the instrument's PSF on the recovered APS is illustrated in the top-left 
panel of Fig. 2: The PSF suppresses the signal at large I and drives it toward the level 
of the Poissonian noise. The typical multipole at which the downturn happens is related 
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Figure 2: Top left: Effect of the PSF on the recovered APS. The slope of the input 
spectrum is s = 0.5. The PSF widths are (Tpsf = 0.1° (red band) and cjpsf = 0.2° (blue 
band). Top right: Effect of the fov on the APS for a pure background (isotropic) event 
list, with cTpsf = 0.05°. The width cTfov is increased in steps of 0.5, in the range from 
Cfov = 0.5° (red) to crfov = 3.0° (blue). Bottom: Influence of the signal fraction /sig 
on the recovered APS for an input slope of s = 0.5; background events are distributed 
isotropically. Dotted lines show the APS in case of a vanishing PSF and fov distortion. 



to the PSF width cjpsf and given by £s w 180°/(Tpsf. This effect can be corrected if cjpsf 
(or the full PSF shape in the non-Gaussian case) is known [1], although at the expense 
of an increasing uncertainty of the recovered APS. In the forward- folding approach that 
we are using here this is, however, not necessary, since the model is directly convolved 
with the PSF before comparison with the (simulated) data (see Section 4). 

Due to a finite fov, anisotropies at a scale larger than the fov will be suppressed. 
This is illustrated in the top-right panel of Fig. 2, where the APS of an isotropic 
(background-only) event list is shown for different fov. A larger fov allows to explore 
larger scales and thus lower multipoles. The minimum resolvable multipole is approxi- 
mately given by £min ~ 180°/(Tfov. 

The bottom panel of Fig. 2, instead, illustrates the effect of different signal-to- 
background fractions /sig. The background is produced mainly by two different processes: 
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(i) Events caused by mostly isotropic cosmic rays (protons and electrons), which have 
been misclassified as photons. Anisotropies in the electron background have not 
been detected so far [35], while anisotropies in the hadronic background are indeed 
present at the level of 10~^ [36-40]. However, they extend only to a multipole of 
~20 (see [36]) and are thus negligible in our analysis. For Cherenkov telescope sys- 
tems the hadron background component dominates, and, therefore, a sufficiently 
good gamma-hadron separation is a crucial characteristic of the instrument. To 
reduce the residual background of cosmic-ray electrons important below a few hun- 
dred GeV, a sufficiently good gamma-electron separation would also be favorable. 
However, gamma-electron separation capabilities are limited, due to the similarity 
of electron and photon initiated showers. With current instruments a rejection of 
~50% of the electrons seems possible [31, 32], while the expected performance of 
CTA on this aspect have not been studied in detail yet. 

(ii) An intrinsically isotropic component of the diffuse photon background, which does 
not act as signal according to our definition. For instance, cosmological photons 
produced by truly diffuse processes may account for this background component. 

Note that we do not consider the effect of the local variations of the night-sky 
background and its effect on the acceptance of the fov [33], assuming that prominent 
features in the fov such as bright stars can be eventually masked and excluded from the 
analysis. A more realistic MC simulation would be required to study the effect in detail, 
which is left for future work. 

In the bottom panel of Fig. 2, a slope of s = 0.5 has been assumed for the input 
spectrum used to simulate the maps and event lists. Dotted lines show the recovered 
APS in the case of vanishing PSF effects as well as fov distortion (but with background 
still included). Clearly, the background-fraction has a large influence and considerably 
reduces the signal height with respect to the intrinsic Poissonian noise level, i.e. the 
signal-to-noise ratio. We conclude that a good background rejection is crucial. 

Finally, we consider the effect of the observation strategy. In general, a large amount 
of observation time is required to detect anisotropies, and a dedicated observation could 
be, eventually, not feasible due to different targets competing for the limited observation 
time available. The use of existing observations obtained on different targets, combining 
them to search for anisotropies, would thus be preferable. This approach is investigated 
in Fig. 3. Here, 10^, 10^, and 10^ events are distributed within one fov and in 10 different 
ones. Qualitatively, when the same number of events is distributed in a single or multiple 
fov two effects are competing. On one hand the Poissonian noise {=Ai}/N in the case 
of partial sky-coverage, where Ail is the covered solid angle) increases linearly with the 
number of fov, decreasing the sensitivity. On the other hand, with multiple fov more 
modes are available for the APS calculation reducing the error on the APS itself (on the 
sphere this effect is known as cosmic variance) and thus increasing the sensitivity. The 
overall number of events determines which error dominates (the Poissonian one or the 
cosmic variance ) and thus which effect is dominant. 
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Figure 3: Influence of the observation strategy on the APS. Shown are the recovered 
APS for an input spectrum with slope s = 0.5 and for 10^ (red), 10^ (green), and 10^ 
(blue) events distributed in a single fov (dark colors) or 10 different fov (light colors). 
In all cases, the APS are normalized to the corresponding noise level, which is shown by 
the solid black lines. The dotted lines show the input APS for the 10^ events cases. 



Fig. 3 illustrates that in the limit of high statistics (the 10 events case), having 
multiple fov indeed decreases the error on the APS, indicating that the measurement 
is cosmic variance limited, and the sensitivity is thus improved. In the limit of low 
statistics (10^ events), instead, the error is approximately unchanged indicating that 
the Poissonian noise is comparable or dominates over the cosmic variance. In this case 
the sensitivity is worsened or remains unchanged. Finally, the 10^ events case is inter- 
mediate giving a marginal improvement in sensitivity. A quantitative numerical study 
of this effect in a realistic scenario with background is presented in the next section. 
Approximate analytical formulae are given in Appendix A. 



3 Benchmark instrumental setups and cosmic-ray backgrounds 

In the following, we simulate two instrument types motivated in the first case by the 
characteristics of the current generation instruments such as H.E.S.S., MAGIC, and 
VERITAS, and in the second case by the expected properties of CTA. We also consider 
two threshold observation energies (100 GeV and 300 GeV) which have different signal- 
to-background ratios (see below and next section). 

For the current generation instruments the performances below 1 TeV typically 
degrade rapidly in energy and we thus consider only a threshold of 300 GeV. Above 
300 GeV, an effective area of 10^ m^ (after selection cuts which improve the fraction of 
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Table 1: Characteristics (effective area Aeff, field of view Cfov, and PSF Cpsf) of the 
instrumental setups that we use as benchmarks to perform simulations. The benchmark 
setups have been chosen in accordance with characteristics of the current generation 
of lACTs, and the planned CTA observatory. Predictions are made for two different 
threshold energies ^Jthr, 100 GeV and 300 GeV. 

gamma rays with respect to hadrons) is typical, together with an angular resolution of 
0.1° and a fov of 1.7° (see, e.g., [23-25, 41]). 

For CTA, recent Monte-Carlo studies of the expected performance [42] (see also 
[26, 43]) indicate an effective area of 3 X 10^ m^ above 300 GeV and 10^ m^ above 100 GeV 
(see figure 15 in [42]). For both cases we assume an angular resolution of 0.05° (see figures 
10 and 17 in [42]). Both effective area and angular resolution are varying functions of the 
energy, but to simplify the simulations we assume them to be constant. The values chosen 
above are the ones close to threshold energies (rather than averages in the energy range 
above the threshold) and thus represent a conservative choice, since the performances 
of the instrument improve with energy. Finally, for cTfov we consider different choices 
between 3° and 5° (see table 3 in [42]). Note that the rather large fov of 0(5°) might 
be provided by the types C, D, and J of the suggested CTA arrays [42]. Moreover, 
an effectively larger fov can be achieved by a dedicated pointing pattern, adjusting the 
pointings of individual telescopes to correspondingly different offsets from the actual 
source position (see figure 3c in [26]). Table 1 summarizes the characteristics of the 
instruments simulated in this study. The performance of H.E.S.S.-II will be in between 
the two setups considered [41]. 

The isotropic hadronic component of the background is dependent on analysis cuts 
and the quality of the gamma-hadron separation. With respect to the present instru- 
ments, CTA will provide some improvement in the performance of the hadron rejection 
process. However, a substantial part of the background is made by cosmic-ray electrons 
which are more difficult to separate from gamma rays. The study in [42] provides a 
simulation of the total expected background from the hadronic and leptonic component, 
that we use to estimate the background in our setups. In particular, we refer to their 
Figure 16 of the integrated background rate per beam above a given energy threshold 
(and for the assumed effective area). 

For a threshold of 100 GeV, the background-rate is in the range 0.01-0.03 Hz, which 
for an angular resolution of ~0.1° and a fov of 5° gives a total background rate in the 
range 50-150 Hz. We will thus use benchmark background rates of 10 Hz and 100 Hz, the 
first being slightly optimistic but still possible depending on eventual improvements of 
the background rejection, and the second more conservative. For a different fov the rates 
scale with the factor (cfov,i/o'fov,2)^- For simplicity, we use 10 and 100 Hz as benchmark 
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background rates for all the considered fov, although it should be kept in mind that for 
fovs of 3° and 4° the rates are over-estimated. Similarly, the background per beam for a 
threshold of 300 GeV is in the range 3 x 10~^-3 x 10"'^ Hz, and the angular resolution in 
the range 0.06°-0.1°, which for a fov of 5° gives background rates in the range 1.5-42 Hz. 
Again, we will discuss an optimistic and more conservative background-rate, i.e., IHz 
and 10 Hz, respectively. 

For the case of current lACTs (H.E.S.S., MAGIC, VERITAS), assuming the same 
characteristics as of CTA above 300 GeV and scaling to the smaller fov, the background 
rate would be reduced by a factor of 10. A further reduction in the rate is given by the 
smaller effective area. On the other hand, the background rejection capabilities are more 
limited with respect to the expectation for CTA, thus increasing the background rate. 
We assume the same benchmark background rates above 300 GeV as for CTA, i.e. 1 Hz 
and 10 Hz. These values are indeed in line with typically observed background rates. 

4 Dark matter sensitivity 

In this section, we use a more realistic setup to simulate maps with a given level of 
anisotropy. Here, we only consider anisotropy with a slope of s = 2, i.e., the same slope of 
the isotropic Poissonian noise. With the conventions given in section 2, s = 2 corresponds 
to an APS constant in multipole and can therefore be characterized by a single number, 
i.e. Ci = Cp. This kind of anisotropy spectrum, known as Poisson anisotropy, is 
typically expected from unresolved point sources, and is a good approximation also 
for most of the DM models, see [7, 8, 10-18]. The similarity to the Poissonian noise 
also suggests a straightforward way to simulate this kind of anisotropy. For A'^ equal 
sources distributed all over the sky the Poisson anisotropy of the fluctuation map will 
be Cp = Att/N , and, inverting the process a map with an anisotropy equal to Cp can be 
simulated distributing N = An/Cp equal sources on the sphere. We will use this method 
in the following. In general, unresolved sources are not identical but have a certain 
flux distribution dN/dS. Indeed, anisotropy measurements can be used to recover the 
underlying dN/dS [2, 5]. A simulation with a delta-like dN/dS is however sufficient for 
the DM sensitivity study presented in the following. Once this source map is produced, 
we follow the algorithm described in Section 2 to produce a sequence of background or 
signal events. Differently from Section 2, the source maps are not rescaled, since they 
possess an intrinsic anisotropy normalization (the one given in input) that we want to 
retain. 

The absolute number of background events is simply given by the integrated back- 
ground rates estimated in the previous section over the observation time. For the total 
number of gamma rays we normalize the flux to the extrapolation above 100 GeV and 
300 GeV of the spectrum of the isotropic gamma-ray background (IGRB) measured with 
Fermi-LAT (l){E) = (/)o (^^/lOOMeV)-^-^!, with (/>o = 1.45 x lO"'^ cm'^ s"! sr"! MeV'^ 
[3], which gives the values (t){E > 100 GeV) 6 x 10"^° cm~^ s"^ sr"^ and ct){E > 
300 GeV) « 1.3 x 10^^° cm"^ sr"^ We neglect the effect of attenuation of the IGRB 
spectrum expected by pair production on the extragalactic background light. The ex- 
pected softening in most of the attenuation models is only mild [44-50], and taking it 
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into account would only slightly reduce the total flux above 100 GeV, while above 300 
GeV the attenuation could be more pronounced. We get 10448 and 6 659 events in 
total (number of events integrated over the entire field of view) for a 100 h observation 
with CTA, assuming a fov of 5° and a threshold of 100 GeV and 300 GeV, respectively 
(setups as discussed in table 1). The events are distributed between DM and astrophys- 
ical sources according to the relative contribution to the IGRB. Notice that while the 
number of gamma-ray events is reduced by less than a factor of 2 from a threshold of 
100 GeV to 300 GeV, the number of background events is reduced by a factor of 10, 
due to their steeper energy spectrum and the improved background rejection at higher 
energies. It is thus convenient to consider a threshold of 300 GeV despite the decreasing 
number of events. The numbers of gamma rays and background events for all the setups 
are reported in table 2. The resulting signal fraction is the fraction between the signal 
rate and the sum of signal and background rate. 

The intrinsic anisotropy of astrophysical sources is modeled in accordance to the 
recent measurement of anisotropy performed with Fermi-LAT in the range 1-50 GeV 
[1]. The fluctuation energy spectrum of the measured anisotropy is compatible with a 
constant value of ~ 10~^, while the intensity energy spectrum of anisotropy is compatible 
with a power law with slope ~ 2.4. Combined, these results indicate that the measured 
anisotropy originates from unresolved blazars. This is further supported by the analysis 
performed in [2]. We assume that the result holds above 100 GeV and we thus use the 
value Cp = 10~^ for the intrinsic astrophysical source anisotropics. We will however 
explore also other values of Cp to assess the robustness of the results with respect to 
the choice of this parameter. 

The theoretical predictions for the intrinsic DM anisotropy are uncertain and span 
over various orders of magnitude from 10~^ to 10"^. We assume a benchmark value 
(jDM _ ]^Q-3 comment on other values. A simplified analytic calculation is 

reported in appendixes A, B and C, which illustrates the expected dependence of the 
sensitivity on the choice of Cp^ and Cp. 

To estimate the sensitivity to the DM component, we vary the relative contribu- 
tion of DM to the total IGRB flux assumed above, assigning the remaining flux to the 
astrophysical component. In particular, we consider five values of the DM fraction p, 
i.e., 2%, 5%, 10%, 22%, and 46%, and we interpolate the results between them. We 
perform 20 realizations, used to estimate the average astrophysical spectrum Cf^'^'^ , its 
error ac^ , and the average spectrum c^^'P'^" composed of both the astrophysical and a 
fractional DM contribution of p. We emphasize that Cf^^"^ ^ 10~^ since the presence of 
the background alters the absolute normalization of the APS. The same consideration 



To quantify the sensitivity to the relative DM contribution, we use a simple chi- 
square approach, comparing two different definitions: 



applies to C^ 



rDM,p% 



^ 10-3. 




(4.1) 
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1000 /^DM,p%,i ^astr 

where C'^'^'^^'-' and C'^^^^'^ are the spectra from the single reaUzations as opposed to the 
average ones, c^^'P'^" and C^'^*'', respectively. In the estimate we use only multipoles 
above 100, discarding lower multipoles affected by the fov of the instrument. The chi- 
square x^(p) follows a ^^"distribution with 901 degrees of freedom, so that we can 
quote sensitivities at 95% confidence level (CL) for the value of p, for which x^{p) 
is larger than 972. However, since scatters comparing different realizations, a 

further criterion must be specified. For example, the average value of xf (p), i.e. Xiip)^ 
can be used. Instead, we use the more conservative requirement that at least 19 out 
of the 20 different xfip) values (for a given p) are larger than the sensitivity threshold 
of 972, in order to set the value of p corresponding to the 95% CL. The two different 
ways to define x^(p)) see Eqs. (4.1) and (4.2), correspond to two different ways to define 
the sensitivity. In the first case, we assume and simulate a "true" case which has no 
DM contribution and we search for the minimum (at 95% CL) DM fraction required to 
exclude the null hypothesis. In the second case, we assume and simulate "true" cases 
with a DM contribution of p and we search for the minimum contribution for which the 
null hypothesis is incompatible with the simulated data. We find the two cases to give 
consistent results. Let us also stress that, in both cases, our definition of sensitivity is in 
short a "95% CL incompatibility with the null hypothesis" , which implies a comparison 
of two x^-distributions. A further commonly employed definition of a sensitivity as "95% 
upper limit in the case of a null detection outcome of the experiment" (which requires 
interval estimation through a profile likelihood or test statistic procedure) can give more 
stringent sensitivities, but will not be used here. 

Finally, we note here that the sensitivity estimation described above assumes a 
previous knowledge of the intrinsic astrophysical and DM anisotropics, while, in realistic 
cases, it would be preferable to measure these quantities from the analysis itself in order 
to perform a comparison with the prediction of the models. This can be done, in practice, 
analyzing the data in multiple bands of energies. The presence of more components at 
different energies and the intrinsic anisotropics of the components can then be inferred 
by a study of the anisotropy energy spectrum in intensity and fiuctuation [13, 17]. We 
do not pursue such detailed analysis here, while for the sensitivity estimate to the DM 
contribution the above approach is sufficiently accurate. 

We show the results in table 2. It is quite clear that current instruments have 
very limited prospects of constraining DM through anisotropics. It is worth stress- 
ing, nonetheless, that looking at anisotropics has also implications for the astrophysical 
sources, as we discuss in more detail below, so that this particular result should not 
hamper a dedicated search for anisotropics with current observatories. However, the 
prospects for DM searches improve for CTA, especially if the background rate can be 
kept reasonably under control, and sensitivities of ~ 15% and ~ 8% at 100 GeV and 
300 GeV, respectively, can be achieved. Even in the most conservative case of a fov of 
3°, still a sensitivity of ~ 30% above 300 GeV can be reached (not shown in the table). 




- 12 - 



H.E.S.S./MAGIC/VERITAS Eth=300 GeV afoy = 1.7' 





Observation time [h] Bg. 


rate [Hz] 


Sens. 




A^bg 








100 




1 

10 




> 46% 

> 46% 


257 


3.6 X lO'^ 
3.6 X 10*^ 








300 




1 

10 




> 46% 

> 46% 


770 


1.08 X 10'' 
1.08 X 10^ 








1000 




1 

10 




> 46% 

> 46% 


2567 


3.6 X 10'= 
3.6 X lO'^ 








10 X 100 




1 

10 




> 46% 

> 46% 


2567 


3.6 X 10" 
3.6 X 10'' 












CTA 


Eth=100 GeV 










Observation time [h] 


Bg. 


rate [Hz] 


ffov ' 

Sens. 




O"fov 

Sens. 


= 5° 

^sig 




TVbg 


100 




10 
100 




>46% 
>46% 


6687 


> 46% 

> 46% 


10448 


3.6 X 10'' 
3.6 X 10' 


300 




10 
100 




30% 
> 46% 


20059 


30% 
> 46% 


31343 


1.08 X 10' 
1.08 X 10® 


1000 




10 
100 




30% 
>46% 


66867 


15% 
> 46% 


104476 


3.6 X 10' 
3.6 X 10* 


10 X 100 




10 
100 




> 46% 
>46% 


66867 


30% 
> 46% 


104476 


3.6 X 10' 
3.6 X 10* 








CTA 


Eth=300 GeV 










Observation time [h] 


Bg. 


rate [Hz] 


ffov 

Sens. 


= 4° 


O'tov 

Sens. 


= 5° 




TVbg 


100 




1 

10 




30% 
>46% 


4262 


30% 
> 46% 


6659 


3.6 X lO'' 
3.6 X 10" 


300 




1 

10 




30% 
>46% 


12785 


15% 
> 46% 


19976 


1 
1 


08 X 10" 
08 X 10' 


1000 




1 

10 




15% 
30% 


42617 


8% 
30% 


66587 


3.6 X 10" 
3.6 X 10' 


10 X 100 




1 

10 




30% 
>46% 


42617 


15% 
> 46% 


66587 


3.6 X 10" 
3.6 X 10' 



Table 2: Sensitivity to detect a self- annihilating DM contribution to the IGRB, utilizing 
anisotropy measurements of its angular power spectrum. The sensitivity is given in terms 
of the minimum detectable DM gamma-ray flux, expressed in percentages of the IGRB. 
The three tables list the sensitivities for instrumental setups similar to current lACTs and 
CTA, respectively, for several observation times, background rates, fov, and observational 
strategies. For reference, the number of simulated signal (A'^^ig) and background (Abg) 
events is listed. A value of > 46% in the table means that the sensitivity is outside our 
tested range. 



while a fov of 4° would give a ~ 15% sensitivity. 

A very interesting result is given by the sensitivity achieved with different ob- 
servation strategies. As shown in table 2, the sensitivity of combining ten different 
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(a) Single fov, lOOOh, Sthr = 100 GeV (b) Multiple fov, 10 x 100 h, Sthr = 100 GcV 




(c) Single fov, 1 000 h, Ethr = 300 GcV (d) Multiple fov, 10 x 100 h, Ethr = 300 GcV 



Figure 4: Comparison between the recovered APS for a pure astropliysical case witli 
Ci = Cp = 10~^ (blue bands) and a case with 46% of the total radiation originating 
from self-annihilating DM with Ce = Cp^ = 10~^ (red bands). An observation with a 
CTA-like telescope system of 1000 h on a single target (left column) and of 10 x 100 h 
splitted on 10 different targets (right column) is considered. The upper row refers to an 
energy threshold of 100 GeV an the lower row to 300 GeV. The two cases in each panel 
refer to background rates of 10 Hz and 100 Hz for 100 GeV, and 1 Hz and 10 Hz for 300 
GeV. The size of the fov is cifov = 5°. The lines describe the estimated noise levels. 



observations of 100 h each is just a factor of 2 worse compared to a continuous 1000 h 
single target observation, and roughly equivalent to a full 300 h observation of a single 
target. Following the discussion at the end of section 2.2, this means that, even with 
CTA we will not be in the regime of statistics where multiple fov produce a reduction 
of the errors, but in the regime where the errors are comparable or slightly larger than 
the single fov observation. The implications of this result are, nonetheless, extremely 
important, since in practice the observation of a single target for 1000 h is practically 
not feasible (apart perhaps the Galactic Center over several years), and the observation 
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of ten different targets for 100 h each is more likely to be realized. Most importantly, 
these observations do not need to be specifically targeted for anisotropy, but observa- 
tions taken with different purposes can be analyzed instead. In this manner, the loss in 
sensitivity of a factor of ~ 2 would be acceptable. 

The other crucial parameters determining the sensitivity are C^^ and C^. To test 
the dependence on these parameters, we performed a further simulation using Cp = 10~'^ 
(instead of Cp = 10~^), keeping the value Cp^'^ = 10~^, and we found that the sensi- 
tivities decrease by a factor ~ 3. This seems to be in good agreement with the analytic 
scaling (Cp/Cp^)^/'^ found in the appendix. Given the strong dependence on these two 
parameters, a firmer prediction on the sensitivity requires pinning down their uncertain- 
ties. More accurate calculations of Cp*^ have been recently reported in [18], indicating 
that Cp^ can be as high as 10~^, dominated by the contribution of the Galactic sub- 
structures over extragalactic ones (see in particular their Fig. 7). Although such a large 
DM anisotropy would push the sensitivity to values better than 1%, the intrinsic emis- 
sion from these very anisotropic structures is expected to be very low, as witnessed 
by the fact that the intensity anisotropics are instead dominated by the extragalactic 
component (see again [18]). A more accurate estimate of Cp, instead, awaits a direct 
measurement with Cherenkov telescopes or further work on the modeling of blazar pop- 
ulations at ~TeV energy. In this respect, the example depicted in figure 4 shows that 
observation of an anisotropy Cp = 10~^ is close to the sensitivity attainable with 1000 h 
of observation time with CTA (for the optimistic background estimate), while it would 
be more challenging for the splitted observation strategy. However, Cp is not precisely 
known, and the possibility of a larger power such as 10~^ would obviously improve its 
detection capabilities, although it would imply a lower DM sensitivity. We thus pro- 
pose explicit observations of these kind of anisotropics even for the current generation of 
instruments. Besides the possibility to constrain a DM contribution, an observation of 
anisotropy would provide a complementary and powerful tool to investigate astrophysical 
TeV sources. 

Finally, we translate the sensitivity given in Table 2 as fraction of the IGRB in 
terms of a more common quantity, the DM self-annihilation cross section, although this 
introduces a further model dependence on the DM signal production model. To normal- 
ize the DM signal we use the cosmological DM model of [51] in its (optimistic) version, 
where DM halos have an NFW profile and the presence of DM subhalos is included (see 
also, e.g., [52, 53]), further boosting the annihilation signal. A less optimistic version 
with no substructures would give a factor of ~ 10 less DM annihilation flux and thus a 10 
times worse sensitivity. On the other hand, using the results of [15] from the Millenium-II 
simulation, an order of magnitude enhancement with respect to the above "optimistic" 
case can be achieved with a correspondingly better sensitivity (see also [54], in particular 
Fig.l). The results can be seen in Fig. 5 for various annihilation channels as function 
of the DM particle mass my.. The figure shows that a sensitivity to the thermal self- 
annihilation cross section {av) = 3 x 10~^^ cm^s~^ can be reached for DM masses up to 
a few hundred GeV. Interestingly, even using the conservative version of the IGRB DM 
model above, sensitivities better than the ones achievable with dwarf spheroidal galaxies 
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Figure 5: Left: Sensitivity of CTA on the DM self-annihilation cross section (av) as 
a function of the DM particle mass m^, for an energy threshold of 100 GeV, 300 h (or 
10 X 100 h) observation time, Cfov = 5°, and 10 Hz background rate (corresponding to 
~ 30% of the IGRB in Table 2). Solid lines correspond to the model of [51] including 
DM annihilation in DM sub-halos, while dashed lines refer to the more conservative 
case of no sub-halos giving a sensitivity approximately a factor of 10 worse. We consider 
different annihilation channels, i.e., bb (red), r^r~ (green), and /i"*"//" (blue) final states. 
See text for more details and for the assumed instrumental characteristics. Right: Same 
as above for an energy threshold of 300 GeV, 300 h (or 10 x 100 h) observation time, and 
1 Hz background rate (corresponding to ~10% of the IGRB in Table 2). 



observed with CTA and comparable with achievable limits using clusters observations 
with CTA can be reached [27]. 

5 Discussion and Conclusions 

With a simplified MC setup, we have investigated the key aspects of the capability 
of ground-based gamma-ray observatories with small fov (i.e. imaging air Cherenkov 
telescopes) to observe anisotropics in the diffuse gamma-ray background. In particular, 
we have studied the effects of the angular resolution (PSF), the field of view, and of the 
hadron-background rejection efficiency, which are found to be crucial characteristics for 
a good sensitivity to detect angular anisotropics. 

Using benchmark detector configurations representing realistic setups for the present 
generation experiments like H.E.S.S., MAGIC, and VERITAS, and the forthcoming 
CTA, and using realistic expectations for the anisotropics from self-annihilating dark 
matter and from astrophysical sources, we have studied the sensitivity to detect a DM 
component through an analysis of angular anisotropics. We find that with an observa- 
tion time of ~1000h and background rates of ~ 10 Hz above 100 GeV and 1 Hz above 
300 GeV, CTA will reach sensitivities of the order of ~ 10% to a relative contribution 
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of DM annihilation to the total isotropic gamma-ray background. More important, we 
find that with a multiple field-of-view observation strategy of 1000 h splitted into ten 
separate targets of 100 h observations yields only a slight reduction in sensitivity (~20%, 
a factor of 2). This means in practice that it will be possible to obtain interesting con- 
strains on DM without specifically dedicated deep observations, but combining instead 
existing observations of different primary astrophysics targets. The sensitivity achievable 
can be already sufficient to test the thermal annihilation cross section for DM masses 
< 200 GeV for certain models of DM annihilation in Galactic and extragalactic environ- 
ments. We also find that CTA will have sufficient sensitivity to detect anisotropics from 
astrophysical sources with a Poisson anisotropy level of 10~^, while the sensitivity of cur- 
rent generation instruments is approximately one order of magnitude lower. Given the 
uncertainty on the exact anisotropy level expected, we nonetheless propose that avail- 
able deep exposures, preferably at high galactic latitudes, should be indeed analyzed to 
search for anisotropics. An observation of anisotropy would provide a complementary 
and invaluable tool to investigate the nature of TeV sources. 
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A Error on the Poisson Anisotropy 

The error on the fluctuation angular power spectrum is given by 



where Cn = ^iov{^/N^ + Nb/N^) is the (Poisson) noise, with A'^^ the number of gamma- 
ray events, Ni, the number of background events, and Qfov the total Field of View (fov) 
in steradians [7]; we = exp((Tpgf £^/2) describes the correction for a Gaussian PSF of 
width (Tpsf (in radians). In the following derivation of the sensitivity, we will assume 
a Poisson-like spectrum, i.e. Ci = Cp. In practice, Cp is estimated by calculating 
the weighted average of the measured angular power spectrum, so that Cp = YlitP(-^i-> 
where = {^/ erf) / l/c^) (the weight of higher multipoles is larger, owing to smaller 
errors, until the PSF error starts to dominate). The corresponding error on Cp is the 
given by {5Cpf = Y^iPli^C^f = 1/(E, V^l). 

A.l Simple calculation 

In the following, we consider both a sufficiently narrow PSF, so that the PSF correction 
is not crucial for the multipole range of interest, as well as the high statistics limit, so 
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that Cm is negligible with respect to Cp. In this limit, we have (implicitly, £ ^ 1 is 
assumed) 



/sky 



(A.2) 



Evaluating the sum in {SCp)"^ = 1/^1) (this can be done analytically or approxi- 

mating the sum as an integral) yields 



5Cp « Cp. 



/sky (^n 



max ^min/ 



Cp 



/sky 



(A.3) 



A.2 More accurate calculation 

Assuming ^ ^ 1, we have 



a I = {Ci + CNwj)^ 



/sky ^ 



(A.4) 



Evaluating the sum {5CpY = 1/(Z^£ l/"^!) (approximating it as an integral) we find 



5Cp 



/sky 



£2 . 

mm 



In ( Cp/Cn + w'^ 



P ■ 



'-^P %sf 



(A.5) 

With £inin ~ 100, cTpsf = 0.05° ^ 1/^min ~ 0.6°, and in the limit of high statistics 
(Cp » Ctv), the expression can be simplified as 



5Cp K, Cp CTpsf 



Cp 



\^ /sky In 



Cp 
eCpf 



/sky 



(A.6) 



which is equivalent to the previous simplified calculation with the definition 
Note that the quantity In ' 



-1 



Cp 



is of order ~ 1 in the high statistics 



limit. In the case discussed in the paper, Cpsf = 0.05° is assumed, so that the first 
condition is satisfied. Further, for a fov of 5° we have Ofov ~ 10~2sr and /sky ~ 10~^. 
For the case of a 1000 h observation and background rate of 10 Hz above 100 GeV we 
have, referring to the number of events in Table 2, Cn « 3 x 10~^, while for 1 Hz above 
300 GeV we have Cat 8 x 10~^, both to be compared with Cp = 10~^. We thus see 
that, even in the most optimistic cases, we are in the condition Cp ~ Cat. Nonetheless, 
for illustration purposes, in the following we will derive results in the regime Cp ^ Cn 
to estimate the sensitivity achievable under optimal conditions. 
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B Sensitivity using the intensity APS 



Calculations are easier in terms of the intensity APS, since power spectra are linearly 
additive in case of uncorrelated summands. This is indeed the case for Poisson-like 
anisotropies. 

The intensity APS is simply related to fluctuation APS as 

Ci, = I^Cp, (B.l) 

where / is the intensity of the component. We consider the same scenario as above, 
i.e., an astrophysical component with Poisson anisotropy Ca and intensity I a, and a 
DM component with anisotropy Cdm and intensity I dm- The corresponding intensity 
anisotropies are given by = I\Ca and C^,^^ = I'jyj^^CoM- Assuming that the 
intensity anisotropy Q^^"^^^^ ± SCp'^^'^^ has been measured in a given energy band, we 
can set a very conservative upper limit on the DM contribution as 

ChM<C'/'''\ (B.2) 

which implies 

T / /^data 



If the DM component is subdominant, we have Cp^^^ ~ Ca and Idm/I ^ Ca/Cdm- 
For our benchmark case with Ca = 10"^ and Com = 10"^ this yields a maximum 
sensitivity of Idm/I ^ 10%. 

A more interesting case, which more accurately approximates the case discussed 
with the simulations above, is when the quantity Cp"^^^^ = Cj^ is known within a certain 
error in advance, either because we have measured this quantity in a different energy 
band, or because it matches our theoretical expectations. In this case, a more interesting 
upper limit can be derived from 



which implies 



ata 



I '"V ^DAI V /sky ^max 



Again, with Cp^*^ = Ca and for our benchmark case Ca = 10~^, Cdm = 10~^, ^max ~ 
1000, and /sky ~ 10~^, we have a maximum sensitivity of Idm/I ^ 2%. Note that the 
sensitivity improves with a larger fov, although only with the fourth root of it. However, 
the sensitivity improves faster with the angular resolution (-^max)* 

We emphasize that the above results are clearly back-of-the-envelope calculations, 
with the purpose of estimating approximately the sensitivity and its dependence on 
the relevant quantities. Dealing with real data, more appropriate derivations can be 
conducted, for example, with a likelihood analysis. 
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C Sensitivity using the fluctuation APS 



In the paper, fluctuation anisotropies are used to calculate the sensitivities. In this case, 
the relation we are effectively imposing to derive the sensitivity is 

^data ^^data <^ (J <^ ^data _j_ ^^data 

where 

[ioM + iA) 

Rewriting 0$^^'" ± (5C|,^*^ as C|,^*'^ (1 ± A), where A = ^2/(/sky£Lx)> we find that the 
sensitivity is given by a second order equation, solved by 



Cdata (1 ± A) - Cdm 



(C.3) 




We can simplify the above equation in a few relevant cases. For example, if Cdm ^ Ca 
(as in the benchmark case), and again assuming Cp^^'"^ ^ Ca, the positive (physical) 
solution is 

/^data / /^data / o 

(C.4) 

For the benchmark numbers we get a sensitivity of Idm/I ^ 3%, only slightly worse 
than in the intensity case. The sensitivity quoted in the paper for the most favorable 
scenario is actually 8%, therefore not very far from the analytical estimate. 

Also note that for Ca = 10""^ and Cdm = 10~^ the above sensitivity degrades quite 
rapidly, getting worse by a factor of 10, in agreement with the simulations performed 
explicitly. In the intensity case, instead, the dependence is only with the square root, 
and the sensitivity should worsen only by a factor of ~3. 

Finally, unlike in the intensity case, we can get an upper limit on I dm, even in the 
case Cdm Ca- This case is quite unphysical, although it can be practically taken 
as case study for a further (non-DM) component with negligible anisotropy. This is 
expected, for example, from normal galaxies or truly diffuse processes like photons from 
UHECRs cascades on the CMB. In this case, we get 



I y /sky^max 

after a few simplifications, which interestingly is independent of Ca- The benchmark 
numbers for /sky and ^max yield a limit of Ida/ /-^ ^ 5%. 

The case Ca ~ Cdm seems to give the worst limit. In this case, a simple formula 
cannot be derived and the full expression Eq. C.3 needs to be used. 
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